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The structural, electronic and magnetic properties of single layers of Iron Dichloride (FeCD) were 
calculated using first principles calculations. We found that the IT phase of the single layer FeCD is 
0.17 eV/unit cell more favorable than its IH phase. The structural stability is confirmed by phonon 
calculations. We found that IT-FeCD possess three Raman-active (130, 179 and 237 cm“^) and 
one Infrared-active (279 cm“^) phonon branches. The electronic band dispersion of the IT-FeCD 
is calculated using both GGA-PBE and DFT-HSE06 functionals. Both functionals reveal that the 
IT-FeGD has a half-metallic ground state with a Curie temperature of 17 K. 

PACS numbers: 81.05.ue, 85.12.de, 68.47.Fg, 68.43.Bc, 68.43.Fg 


The synthesis of stable two-dimensional (2D) 
materiala^^ and observing their various exceptional 
physical properties triggered a large interest in their 
potential technological applications. Single layer 
transition-metal dichalcogenides (TMDs) are a new 
class of materials which are promising candidates for 
next generation of flexible nanoelectronic devices due to 
their wide range of important properties in their bulk 
and monolayer form such as superconductivity^ and 
half-metallicity<^“— 

The single layer TMDs can be obtained from their 
three-dimensional (3D) counterpart which consist of 
weakly interacting 2D layered structures. Most of these 
materials have either De/ (2H structure) or Ds^ (IT 
structure) point-group symmetry and very few of them 
are stable in both the IT and 2H phases. A lot of efforts 
are being done to predict new 2D materials theoretically 
and synthesize them experimentally with tailored elec¬ 
tronic, optical, and chemical properties. 

Half-metallicity is one of the desired properties for fu¬ 
ture spintronics devices and it is observed in 3D com¬ 
pounds several decades ago.— Most of the 2D materi¬ 
als do not have intrinsic half-metallicity but it has been 
reported that graphen e^ and hydrogenated boron- 
nitride nanoribbonsH gain half-metallicity when in-plane 
electric fields were applied. Chemical modification 
of the edges with different functional groups is also 
an alternative approach to achieve half-metallicity in 
2D structures . It has been recently reported that 
graphitic carbon nitride (g-C 4 N 3 ) ^^d'^ , Mn 02 nanosheet 
with vacancies and monolayer chromium nitride^^ show 
half-metallic behaviour. 

FeCb, which is a metal in its bulk form, is another 
member of family of TM crystals. Its bulk and molecular 
form have been investigated. Saraireh and Altarawneh 
— investigated the structure and thermodynamics sta¬ 
bility of the FeCl 2 surface using DFT calculations and 
found that the FeCl 2 (100-Cl) surface is the most sta¬ 
ble configuration. The electronic energy levels and the 
Mossbauer spectroscopic parameters of FeCl 2 were cal¬ 
culated by Bominaar et alj^ and by Mishra et ali^ at 
the self-consistent field level. The X-ray properties of 
FeCl 2 were reported by Chou et al— and Veal et al— us¬ 



FIG. 1: (Color online) Atomic structure of 2D single-layer 
FeGb in the T structure, (a) Top and (b) side views of the 
T structure showing the primitive unit cell of the 2D hexag¬ 
onal lattice with Bravais lattice vectors a and b (|a| = |b| and 
relevant internal structural parameters. 

ing density functional theory. Shechter et al. performed 
a Mossbauer resonance study on FeCl 2 monolayers de¬ 
posited on oriented basal planes of graphite— 

In this work we will investigate the structural, elec¬ 
tronic and vibrational properties of monolayer FeCl 2 by 
using first principle calculations. We will show that the 
minimum energy configuration of the single layer FeCl 2 
is the T structure which has a half-metallic behaviour. 

All the calculations are performed using the projec¬ 
tor augmented wave (PAW)^ potentials as implemented 
in the Vienna Ab-initio Simulation Package (VASP)— >22 
The electronic exchange-correlation potential is treated 
within the spin polarized generalized gradient approxi¬ 
mation (GGA) of Perdew-Burke-Ernzerhof (PBE)— . A 
plane-wave basis set with kinetic energy cutoff of 500 eV 
is used. A vacuum spacing of more than 15 A is taken to 
prevent layer-layer interactions. A set of (20x20x1) F 
centered k-point sampling is used for the primitive unit 
cell and scaled according to the sizes of the supercells 
used. The structures are relaxed until self-consistency 
is reached at the level of 10“® eV between two consecu¬ 
tive steps and the atomic Hellman-Feynman force is less 
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than 10“"^ eV/A. Pressures on the lattice unit cell are de¬ 
creased to values less than 1.0 kBar. For the charge trans¬ 
fer analysis, the effective charge on atoms is obtained by 
the Bader methodi^ More accurate electronic structure 
calculations were performed using the screened-nonlocal- 
exchange Heyd-Scuseria-Ernzerhof (HSE) functional of 
the generalized Kohn-Sham scheme^^*^ 

Experimental and theoretical studies have shown that 
most of the compounds of transition metals form sin¬ 
gle layer crystal structures in T (TaS 2 , VS 2 HfSe 2 ) or H 
(M 0 S 2 , WSe 2 , MoTe 2 ) phase. Our calculations show that 
the minimum energy configuration of a single layer FeCl 2 
is the IT structure. IH structure of FeCl 2 is 0.17 eV/unit 
cell higher in energy than the IT structure. In the IT 
configuration, the Fe atoms are sandwiched between lay¬ 
ers of trigonally arranged Cl atoms and each Fe atom 
is coordinated to six Cl atoms. The optimized atomic 
structure of hexagonal FeCl 2 lattice is shown in Fig. [1] 
and corresponding structural parameters are presented 
in TableHl The calculated lattice parameter, 3.47 A, is in 
good agreement with previously reported experimental 
and theoretical values which are around 3.5 AA^i^ How¬ 
ever, the optimized lattice parameter of the H phase is 
3.34 A which is smaller than the one for the T phase. The 
Fe-Cl bond length is 2.45 A in the T phase and is almost 
the same in the H phase which is 2.47 A. Another possi¬ 
ble phase of monolayer TMDs is the distorted T structure 
which was observed in ReS 2 i^ In ReS 2 , expanding the 
lattice parameter of the IT-ReS 2 in one direction leads to 
a dimerization of Re atoms. This is responsible for an en¬ 
ergetically more favorable distorted T structure. A sim¬ 
ilar investigation for the FeCl 2 monolayer revealed that 
the distorted T phase is not likely in this compound. This 
observation fits with earlier theoretical and experimental 
works 30,32 

The cohesive energy per unit cell of FeCl 2 is calculated 
by subtracting the total energies of the individual atoms 
from the total energy of FeCl 2 which is formulated as 
Ec = EtIFc] + 2Et[CI] — ET[FeCl 2 ]- We found that the 
cohesive energy of FeCl 2 per unit cell is 2.64 eV/atom in 
the T phase which indicates a strong cohesion relative to 
the constituent free atoms and it is slightly higher than 
the value for the H structure which is 2.59 eV/atom. 

TABLE I: Calculated values of stable, free-standing, 2D 
single-layer FeCb in the T structure: lattice constant, |a| = |b|; 
bond lengths, dpe-ci and dci-cv, Cl-Fe-Cl bond angle, 6\ to¬ 
tal magnetic moment per the unit cell, p; donated (receieved) 
electron by per Fe (Cl) atoms, ppe (pci)- 


FeCL a 

dpe-Cl 

dci-ci 

e 


pPs 

Pci 

E, 

Phase (A) 

(A) 

(A) 

(deg) (/is) 

(e) 

(e) 

eV/atom 

IT 3.47 

2.45 

2.84 

90.15 

4 

-1.20 

0.60 

2.64 

IH 3.34 

2.47 

3.08 

77.23 

4 

-0.86 

0.43 

2.59 


FeCl 2 is not the only possible FeCl„ type structure in 
the monolayer form. Indeed, recently it was reported 
that FeCla also exists<^ In order to investigate the most 
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FIG. 2: (Color online) Phonon dispersion for the IT phase of 
FeCb. Eigenfrequencies for Raman and Infrared-active modes 
and corresponding eigenmotions are shown in the lower panel. 

favorable chemical composition of the compound, we cal¬ 
culated the cohesive energy of FeCls. The calculated co¬ 
hesive energy of FeCls is 2.4 eV/atom which is smaller 
than the value for the IH and IT phases of FeCl 2 . This 
shows that the IT phase of FeCl 2 is the most favorable 
monolayer structure of the FeCl„ type structures. 

Phonon dispersion of the single layer FeCl 2 in the IT 
phase is shown in Fig. [21 Here the dynamical matrix and 
the vibrational modes were calculated using the small- 
displacement method (SDM)^ with forces obtained from 
VASP. We found that all the phonon branches have real 
eigenfrequencies through all the symmetry points and 
hence the predicted structure of lT-FeCl 2 is dynamically 
stable. 

There is one Fe atom and two Cl atoms within the 
primitive unit cell of IT-FeCl 2 and its phonon spec¬ 
trum includes 9 phonon bands, 3 acoustic (transverse 
acoustic TA, longitudinal acoustic LA, and out-of-plane 
transverse acoustic ZA) and 6 optical. While LA and 
TA acoustic branches have characteristics linear disper¬ 
sion, the frequency of the out of plane ZA mode has a 
quadratic dispersion in the vicinity of the zone center. 
Analysis of the lattice dynamics shows that the decom¬ 
position of the vibration representation of optical modes 
(without translational acoustic modes) at the F point is 
Topt = 2E' -I- Ai' -I- 2E" -I- A 2 ". We also show the optical 
character and frequency of the phonons for the mono- 
layer FeCl 2 in the lower panel of Fig. |21 Note that while 
the phonon modes that have out-of-plane character {Ai 
and A 2 ") are singly degenerate, modes with in-plane vi¬ 
brational character {E' and E") are doubly degenerate. 

In addition a large phonon energy bandgap, which is 
130 cm“^ at the F and 30 cm“^ at the K point, is 
found between the acoustic phonon modes and the op- 
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FIG. 3: (Color online) Electronic band structnre of 2D single¬ 
layer FeCD in the T structure and projected DOS of Fe and 
Cl atoms (a) minority and (b) majority states. The Fermi 
level of both calculations are set to 0 eV. 


tical phonon modes. Due to the large phonon bandgap, 
the acoustic vibration is preserved from being interrupted 
by the optical modes and hence, FeCl 2 may have a higher 
quality factor than a graphene resonator. 

FeCl 2 has a polar character in addition to the cova¬ 
lency of the bonds. Effective charge on cation and an¬ 
ion, Z* and Z*, charge transferred from cation to anion, 
Sp = Z* — Zy [Zy being the valency of the constituent 
atom) are calculated using the Bader analysis. In spite 
of the ambiguities in finding the true effective charge, 
the calculated effective charges in Table H] give some idea 
about the direction of charge transfer and ionicity of the 
honeycomb structure. Bader charge analysis showed that 
the Fe atom donates 1.20 e and each Cl atom receives 0.60 
e charge in the T structure. In the H structure, however, 
Fe atom donates 0.86 e charge and each Cl atom receives 
0.43 e charge. 

In Fig. |3]the band structure and the density of states 
of FeCl 2 are shown for the minority (Fig. [Ha,) and the 
majority (Fig. (Hd) states. The states which are around 
the Fermi level are mostly from the d-electrons of the Fe 
atom. The bands which are plotted in black in Fig. [Hare 
the results of PBE-GGA calculation. According to PBE- 
GGA calculations EeCl 2 show half-metallic property. An 
energy gap around 4.4 eV is observed in the electronic 


structure of the minority states, however there are bands 
which are crossing the Fermi level in the case of the ma¬ 
jority states which makes the compound a half metal. On 
top of PBE-GGA calculations we also performed HSE06 
hybrid functional calculation to interpret the electronic 
structure of the compound more accurately. The HSE06 
bands are shown in red in the band structure plots. As 
can be seen from the figure the half metallic behavior of 
EeGl 2 remains in the HSE06 calculations and the gap for 
the minority electrons increased to 6.7 eV. 

Ferromagnetically ordered magnetic moments are a 
crucial requirement for spintronics devices. To find 
whether or not ferromagnetic state is energetically fa¬ 
vorable, we perform supercell calculations and compare 
the energies of ferromagnetic and anti-ferromagnetic con¬ 
figurations. For FeGl 2 we compared the total energies of 
non-magnetic, ferromagnetic and anti-ferromagnetic con¬ 
figurations. The spin density plots of these calculations 
are shown in Fig. [H Since the primitive cell of FeGl 2 
has 4 /is magnetic moments, ferromagnetic configura¬ 
tion of 2 X 2 supercell (Fig. |4^) carries 16 pB- The 
anti-ferromagnetic configuration of 2 x 2 supercell (Fig. 
|3 }d) has 0 pb magnetic moment and 93 meV per primi¬ 
tive cell higher energy than the ferromagnetically ordered 
one. The non-magnetic configuration of FeGl 2 is 1.3 eV 
per primitive cell higher in energy than the ferromagnetic 
one. This means that the ferromagnetically ordered state 
is energetically the most favorable arrangement of FeGl 2 . 

Gurie temperature is defined as the temperature which 
separates the paramagnetic phase in which the mag¬ 
netic moments are randomly oriented from the ordered 
ferromagnetic phase. The Gurie temperature of the 
materials can be approximated using the Heisenberg 
model in which the Hamiltonian can be written as H = 
— Yhi j J’mi.ihj, where J is the Heisenberg exchange pa¬ 
rameter and m is the magnetic moment of each site in 
Pb- The expression for the Heisenberg exchange param¬ 
eter for our system is J = (l/12)£'ea;/2m^ where Eex 
is the energy difference between ferromagnetic and anti¬ 
ferromagnetic configurations and the factor (1/12) is due 
to the double counting of the exchange interaction of 6 
nearest neighbor atoms in the summation. If we substi¬ 
tute m = 4 and E^c = 372 into the equation, we get 
J = 0.97 meV. There exist several approximations to 
find Tc from the J value. Here we will use the mean 
field approximation which is kBTc = | J. This approx¬ 
imation leads to a Gurie temperature of 17 K. In ad¬ 
dition, it is worthwhile to note that one may find dif¬ 
ferent, but in the same order, Tc values using different 
methodologiesi^i^i^^. 

To summarize, in this paper we investigated the struc¬ 
tural, electronic and magnetic properties of single layers 
of FeGl 2 using first principles calculations. We found 
that the T phase is energetically more favorable than 
the H phase of the compound. Our vibrational analy¬ 
sis showed that all the phonon branches have real eigen¬ 
values through all symmetry points which indicates the 
stability of the lT-FeGl 2 structure. The electronic struc- 






































4 


a) b) 


FIG. 4: (Color online) Spin density plot of (a) ferromagnetic 
and b) anti-ferromagnetic configurations of FeCb- The blue 
lobes correspond to the minority and the red ones to the ma¬ 
jority spin density. 


ture of lT-FeCl 2 is investigated using both GGA-PBE 
and DFT-HSE06 functionals. A very interesting physi¬ 
cal property of the compound is revealed with both func¬ 
tionals: FeCl 2 has a half-metallic nature. Our super¬ 
cell calculations showed that the ferromagnetic ordering 
of the magnetic moments have a lower energy than the 
anti-ferromagnetic ordered state. Thus, the lT-FeCl 2 is 
an intrinsic half-metallic ferromagnet which makes it an 
attractive compound for future spintronics devices. 
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